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Abstract: This paper proposes a novel, simple and minimally invasive method 
for stroke volume variation assessment using arterial blood pressure measure- 
ments. The arterial blood pressure signal is reconstructed using a semi-classical 
signal analysis method allowing the computation of a parameter, called the first 
systolic invariant INVSi. We show that INVSi is linearly related to stroke 
volume. To validate this approach, a statistical comparaison between INVSi 
and stroke volume measured with the PiCCO technique was performed during a 
15-mn recording in 21 mechanically ventilated patients in intensive care. In 94% 
of the whole recordings, a strong correlation was estimated by cross-correlation 
analysis (mean coe5icient=0.9) and linear regression (mean coefHcient=0.89). 
Once the linear relation had been verified, a Bland-Altman test showed the 
very good agreement between the two approaches and their interchangeability. 
For the remaining 6%, INVSi and the PiCCO stroke volume were not corre- 
lated at all, and this discrepancy, interpreted with the help of mean pressure, 
heart rate and peripheral vascular resistances, was in favor of INVSi. 

Key-words: Arterial blood pressure, first systolic invariant, PiCCO, semi- 
classical signal analysis, stroke volume variation 
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Validation d'une nouvelle methode pour 
I'estimation du volume d'ejection systolique: 
comparaison avec le PiCCO 

Resume : Get article propose une nouvelle methode pour I'estimation du vol- 
ume d'ejection systolique par des mesures de pression arterielle. Le signal do 
pression est reconstruit a I'aide d'une methode d'analyse semi-classique perme- 
ttant le calcul d'un parametre, appele le premier invariant systolique INVSi. 
On montre que INVSi est lincairemcnt relic au volume d'ejection systolique. 
Afin de valider cette approche, une comparaison statistique entre INVSi et 
le volume d'ejection systolique mesure par la technique PiCCO a ete eflFectuee 
pour un enregistrement de 15 minutes pour 21 patients mccaniqucment ventiles 
et en soins intensifs. Pour 94% de I'enregistrement complet, une forte correla- 
tion a ete estimee par une analyse cross-correlation (coefficient moyen=0.9) et 
une regression lineairc (coefficient moycn =0.89). Une fois la relation lineaire 
verifiee, un test de Bland- Altman a montre une bonne correspondance entre les 
deux approches et leur interchangeabilite. Pour les 6% restant, INVSi et le 
volume d'ejection calcule par PiCCO n'ont pas etc correles, et cette difference, 
interpretee a I'aide de la pression moyenne, de la frequence cardiaque et des 
resistances vasculaires peripheriques a ete en faveur de INVSi ■ 

Mots-cles : Pression arterielle, premier invariant systolique, PiCCO, analyse 
semi-classique du signal, variations du volume d'ejection 
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1 Introduction 

Hemodynamic monitoring is crucial for critical care patient management. A 
recent international consensus conference recommended against the routine use 
of static preloaded measurements alone to predict fluid responsiveness [2], and 
dynamic assessment now seems more useful. Several studies have documented 
the ability of respiratory stroke volume variation (SVV) to predict the fluid 
responsiveness in hemodynamically compromised patients [5]- [IS]- Measuring 
respiratory SVV requires a continuous monitoring of stroke volume (SV) which 
can be obtained using invasive or non-invasive methods. Current invasive meth- 
ods used have the disadvantage of requiring the insertion of a central venous 
catheter and the calibration of the cardiac output measure with a cold iso- 
tonic sodium chloride bolus (PiCCO technology) or a lithium chloride bolus 
(LiDCO technology) [9J. An alternative method, which does not require ve- 
nous catheter insertion or calibration has been proposed (Flotrac Vigileo) [5J, 
but several clinical studies have pointed out its poor agreement with reference 
techniques |12| . |16j . Esophageal echo-doppler is the main non- invasive method, 
calculating aortic blood flow from the echo-derived aortic diameter and the 
doppler-derived aortic blood velocity [TT]. Nevertheless, this technique has po- 
tential contraindications, such as esophageal varices or esophageal surgery, and 
several limitations: for instance, it measures the blood flow in the descending 
aorta and not the whole cardiac output. Moreover, the precision of the mea- 
surement depends on accurate probe positioning, which is not always easy to 
obtain [4J. Thus, each of the above methods has its own drawbacks, and there is 
still a need for an easily applicable, minimally invasive, accurate and affordable 
method to estimate SVV. 

Due to the fact that Arterial Blood Pressure (ABP) can be measured us- 
ing minimally invasive or noninvasive methods, the idea of estimating SV from 
ABP has captured scientists for a long time. Thus, many methods have been 
developed and whose objective is to find a relation between one or several pa- 
rameters characterizing the shape of the pressure and SV or cardiac output 
(CO), see for instance [7], [15], [22], [23] and the references quoted there. These 
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methods, which are based on some models of systemic circulation, are called 
pulse contour methods. A comparison between some of the pulse contour meth- 
ods has been proposed in [1 , [26^, ^23, [55]. The simplest model supposes a 
proportionality between CO and the Mean Arterial Pressure {MAP). Other 
approaches, based on windkessel models, link SV to different lumped param- 
eters such as pulse pressure, the systolic and diastolic pressures [7J. However 
these approaches consider the arterial system as a lumped system which appears 
not sufficiently accurate. So, other methods resulting from distributed arterial 
models use the pressure area so that SV is often supposed to be proportional to 
the area under the systolic part of the pressure curve. Corrected versions of this 
relation have been also proposed [15 . However, this approach requires detecting 
the end of the systole which is completely nontrivial, particularly in peripheral 
ABP waveforms. Moreover, approaches taking into account the nonlinear as- 
pects of the arterial system have been proposed, for example modelflow [28], 
but some studies have revealed the poor efficiency of this method in a number 
of cases [24] . 

In this paper we introduce a novel technique for SVV assessment using 
ABP measurements. This method is based on the analysis of ABP with a new 
signal analysis method that was recently proposed in [5T], and called Semi- 
Classical Signal Analysis (SCSA). The new spectral parameters provided by 
SCSA, eigenvalues and invariants, have already given promising results in some 
other applications, as summarized in the following. 

On the one hand, we assessed their ability to discriminate between differ- 
ent situations. In the first situation, nine heart failure subjects were compared 
to nine healthy subjects. In the second situation, eight highly fit triathletes 
were compared before and after training. SCSA parameters always provided 
more significant results than classical parameters, regarding temporal as well 
as spectral parameters ([lOJi [HD- On the other hand, we tested the ability of 
the invariants to represent physiological parameters of great interest, particu- 
larly SVV , in two well-known conditions: the head-up 60 degrees tilt-test and 
the handgrip-test Let us focus on the first invariants. The first global 

invariant [INVi) is, by definition, the mean value of the ABP signal, which is 
a standard parameter in clinical practice. The first systolic {INVSi) and dias- 
tolic {INVDi) invariants are less obvious. They result from the decomposition 
of the pressure into its systolic and diastolic parts. In particular, INVSi corre- 
sponds to the integral of the estimated systolic pressure with SCSA. Referring 
to the pulse contour method stating that the area under the systolic part of 
the pressure curve is proportional to SV as described above, one can show that 
INVSi variations give information on SVV . 

We study in this paper the correlation between INVSi and measured SVV 
using a reference method; the PiCCO technique. The PiCCO technique uses 
the pulse contour method with a calibration by a transpulmonary thermodi- 
lution and is considered a reliable technique. In what follows, we present the 
experimental protocol and recall some basic aspects of the SCSA method. We 
introduce INVSi and its relation to SVV. Then, we present statistical results 
on 21 patients' recordings. 
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2 Materials and Methods 

This prospective study was conducted in the 16-bed medical-surgical intensive 
care unit (ICU) of the Sud-Francilien General Hospital (Evry, France). 

2.1 Patients 

• Inclusion criterion: all mechanically ventilated patients whose cardiac 
output was continuously monitored with a transpulmonary thermodilu- 

tion catheter (PiCCO, Pulsion Medical Systems, Munich, Germany) were 
included, except those satisfying the following excluding criteria. PiCCO 
is routinely used in this unit to monitor hemodynamically compromised 
patients. 

• Exclusion criteria: patients presenting cardiac arrhythmias or breath- 
ing spontaneously were excluded because the SVV is not applicable for 
such patients. 

• Protocol: all patients were sedated with midazolam and fentanyl in 

dosages that wore titrated to achieve full adaptation to the ventilator. 
Ventilator settings were as follows: volume assist-control mode; tidal vol- 
ume (Vt), 6ml /kg ideal body weight; breathing rate, 20 cycles/minute; 

inspiratory /expiratory ratio, -; and Fi02 adjusted to maintain transcu- 
taneous oxygen saturation in blood 94%. Positive end-expiratory pressure 
(PEEP) was set at 5cm H20 but some hypoxemic patients required an 
increase in PEEP to 10 cm H20 during the data acquisition, to improve 
arterial oxygenation. The increase in PEEP was left to the discretion 
of the attending physician, as well as the adaptation of vasoactive drugs 
dosages, adjusted to maintain an adequate circulatory status during the 
protocol. 

2.2 Data acquisition 

One-lead electrocardiogram, arterial pressure, and respiratory flow signals were 
recorded during a 15-min period using a Biopac 100 system (Biopac systems, 
Goleta, CA, USA). All data were sampled at IQQQHz and stored on a hard 
disk. Cardiac output was calibrated just before the data acquisition with a cold 
isotonic sodium chloride bolus of 20 ml. Then, CO and peripheral vascular 
resistances {PVR) were delivered every 30 seconds during the 15-min period. 

2.3 Signal analysis 

Signal processing was performed using the Scilab and Matlab environments at 
the French National institute for Research in Computer Science and Control 
(INRIA-Sisyphe team). 

A Semi-Classical Signal Analysis method 

In this section, wc introduce the SCSA technique and some results of its appli- 
cation to ABP analysis. We also show the relation between INVS\ and SVV. 
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The SCSA principle Let y : 1 1 — > y{t) be a real valued function representing 
the signal to be analyzed such that: 

y G i}(M), y{t) >0,Vte R, 

^eL\Rlm = 1,2, (1) 

with, 

/ + 00 
\Vmi + \t\)dt<oo}. (2) 
-oo 

The main idea in the SCSA consists in interpreting the signal y as a mul- 
tiplication operator, — ^ y.cj), on some function space. Then, instead of the 
standard Fourier Transform, we use the spectrum of a regularized version of 
this operator, known as the Schrodinger operator in for the analysis of 

y- 

H{h;y)^~h'^^y, (3) 

for a small h > 0. The SCSA method is better suited to the analysis of some 
pulse shaped signals than the Fourier Transform |21) . 

In this approach, the signal is a potential of the Schrodinger operator H{h; y). 
We are interested in the spectral problem of this operator which is given by: 

-e'^-y^^X^, telR, (4) 

where A, A G M and V, -0 G H'^{R) Q are respectively the eigenvalues of H{h; y) 
and the associated eigenfunctions. Under equation ([T]), the spectrum of H{h; y) 
consists of: 

• a continuous spectrum A > 0, 

• a discrete spectrum composed of negative eigenvalues. There is a non-zero, 
finite number Nh of negative eigenvalues of the operator H{h;y). We put 
A = with Knh > and ki^ > K2h > • • ■ > Unh, n ^ 1, • • • , Nh- Let 
tpnh, n = 1, - ■ ■ , Nil be the associated L^-normalized eigenfunctions [21] . 



The SCSA technique consists in reconstructing the signal y with the discrete 
spectrum of H(h; y) using the following formula: 

yH{t)^AhY,^^nhi'lh{t), teR. (5) 

Here, the parameter h plays an important role. As h decreases, the approxi- 
mation of the signal improves. However, as h decreases, the number of negative 
eigenvalues Nh increases and hence the time required to perform the compu- 
tation increases. So, in practice, what we are looking for is a value of h that 
provides a sufficiently small estimation error with a reduced number of negative 
eigenvalues. We summarize the main steps for reconstructing a signal with the 
SCSA as follows [H]: 

^_ff^(IR) denotes the Sobolev space of order 2 
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1. Interpret the signal to be analyzed ?/ as a potential of the Schrodinger 
operator H{h; y) ^ ; 

2. compute the negative eigenvalues and the associated L^-normalized eigen- 
functions of H{h; y) ; 

3. compute according to equation ([5]) ; 

4. look for a value of h to obtain a good approximation with a small number 
of negative eigenvalues. 



ABP analysis with the SCSA Now, we introduce some results on the ap- 
plication of the SCSA to ABP analysis. We denote by P the ABP signal and P 
its estimation with the SCSA such that: 

P{t)=4hY,^nh€hit), (6) 
n=l 

where —K.'^h^ n = 1, ■ ■ ■ , Nh are the Nh negative eigenvalues of the Schrodinger 
operator H(h;P) and ipnh the associated L^— normalized eigenfunctions. 

The ABP signal was estimated for several values of the parameter h and 
hence Nh- Fig[T] illustrates measured and estimated pressures for one beat of an 
ABP signal and the estimated error with Nh — 9. Signals measured at the aorta 
(invasively) and at the finger (non invasively) respectively were considered. We 
point out that 5 to 9 negative eigenvalues are sufficient for a good estimation of 
an ABP beat |fT7]. [T^. 

One application of the SCSA to ABP signals consists in decomposing the 
signal into its systolic and diastolic parts. This application was inspired by 
a reduced model of ABP based on solitons solutions of a Korteweg-de Vries 
(KdV) equation d proposed in [TU], [IH|. As described in [T7], [H], the idea 
consists in decomposing ^ into two partial sums: the first one, composed of 
the Ns {Ns = 1,2,3 in general) largest nnh and the second composed of the 
remaining components. Then, the first partial sum represents rapid phenomena 
that predominate during the systolic phase and the second one describes slow 
phenomena of the diastolic phase. We denote by Ps and Pd the systolic pressure 
and the diastolic pressure respectively estimated with the SCSA. Then we have: 

PS) = ^hY,n^hi>lh{t), (7) 

n=l 

Pdit) ^ 4h J2 ^nh^lhit)- (8) 

n=Ns + l 

Fig[2]shows measured pressure and estimated systolic and diastolic pressures 
respectively. We notice that Ps and Pd are respectively localized during the sys- 
tole and the diastole. 



■^Solitons are solutions of some non-linear partial derivative equations like the KdV equation 
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- Estimated pressure 

- Measured pressure 
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Figure 1: Estimation of the pressure at the aorta and the finger level with the 
SCSA and — 9. On the left, the estimated and measured pressures. On the 
right, the relative error 
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Figure 2: (a) Estimated systolic pressure, (b) Estimated diastolic pressure 



INRIA 



Validation of a New Method for Stroke Volume Variation Assessment 



9 



SCSA parameters As seen previously, the SCSA technique provides a new 
description of the ABP signal with some spectral parameters which are the 
negative eigenvalues and the so called invariants [^. The latter consist in some 
momentums of the Unh, n ^ 1, ■ • • , Nh- So we define the first two global invari- 
ants by: 

INVi = 4/1 ^ Knh. INV2 = Y^T. ""nh- (9) 
n— 1 n—1 

Systolic {INVSi,2) and diastolic {INV Di^2) invariants are deduced from 
the decomposition of the pressure into its systolic and diastolic parts and are 
then given by: 

Ns 

INVSi=AhY^Knh, (10) 
n=l 

16 

INVS2^-hY,<h, (11) 

n=l 

INVDi=Ah J2 ""^h, (12) 

n=Ns + l 

16 

INVD2^-h ^ nlf^. (13) 

n=Ns + l 



INVSi for SW estimation We will see here how INVSi is related to 
SV. For this purpose, we recall one approach of the pulse contour methods 
that supposes proportionality between SV and the area under the systolic part 
of the pressure curve as described in the introduction. We denote this area by 
Psa (see figin])- So we have: 

SVpc = kPsa, (14) 

where fc is a positive real and SVpc is the stroke volume estimated with the 
pulse contour method. 

Referring to ([7]) and (|10p . and remembering that %l!nh are L^-normalized, we 
have: 

/ + CX) 
Ps{t)dt. (15) 
-00 

So, INVSi refers to the area under the systolic curve P,. Thus, one can remark 
that both Psa and INVSi describe the area under the systolic pressure but they 
may not be equal because the detection of the end systole in the two cases is 
not the same (see figUja and figE]). Indeed, while Pga is computed by detecting 
the dicrotic notch which is completely non trivial in peripheral ABP waves, 
INVSi results from a nonlinear model of ABP based on solitons that considers 
the propagation of the pulse wave [10], [18] as was described in section Il.C.l.b. 
We can write the following relation between the two areas: 

INVSi - Psa + h, (16) 

^We call these parameters invariants because they are related to the Korteweg-de Vries 
invariants in time 
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Figure 3: Area under the systolic part of the pressure curve used to estimate 
SV 



where 6 G M represents the difference between the two areas. Then, we get: 

INVSi = aSVpc + h. (17) 

with a = Hence, INVSi and SVpc are Unearly related. 
k 

Resulting time series used to compare INVSi to PiCCO stroke vol- 
ume 

On top of INVSi , two vascular time series were analyzed: the heart rate (HR) 
was computed from the pulse interval (PI), which is the distance between two 
systolic occurrences; MAP was calculated from the systolic and diastolic val- 
ues. All data were resampled at 4 Hz, by the interpolation of a third order 
spline function to obtain equidistant data and to guarantee their synchroniza- 
tion. They were then averaged over 15 seconds and delivered every 30 seconds, 
like the PiCCO protocol. PiCCO cardiac output was divided by HR to give a 
stroke volume (SVpicco)- According to the relation between INVSi and SV, 
INVSi was subsequently called SVscsA when it was estimated with the linear 
equation (fT7|) : SVscsA — aSVpicco + i>- Thus, in addition to the SVpicco 
and SVscsA^ temporal relations with HR, MAP and PVR time series were an- 
alyzed to help interpret SV behavior in case of a divergence between SVpicco 
and SVscSA- 

2.4 Statistical analysis 
1. Cross-correlation analysis. 

The cross-correlation analyzes the temporal similarity between two time 
series by estimating the correlation between one time series at time t and 
the other at time t ± x lags (in samples) [5S]- A cross-correlation was 
performed between SVpicco and INVSi time series averaged every 30 
seconds. Cross-correlation coefhcients were computed using the Matlab 
xcorr function (The MathWorks, Inc) after subtracting the means from 
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the time series. Cross-correlation coefficients were computed for all lags 
(—30; +30), each lag corresponding to a 30-second interval. The correla- 
tion coefficients of the unlagged data stand at the midpoint (lag 0). A 
5% level of probability for the correlation coefficients was considered sig- 
nificant (Bravais-Pearson table). An average estimate of the correlation 
over all the subjects was allowed after homogeneity tests (non significant 
z-test and Jarque-Bera test) [T3]. The individual correlation coefficients 
were averaged for each lag. 

2. Linear regression. 

According to the proportional relation between SVpicco and INVSi 
(equation (jl7p '). a linear regression analysis was applied, using SigmaS- 
tat. This analysis provides the Pearson R coefficient, which measures the 
degree of linear correlation between the two estimates, and the parameters 
of the linear equation (jl7l) . a and b. So, they allow an estimation of the 
stroke volume SVscsA using ([T7)) . This transformation is required before 
using the Bland- Altman test. 

3. Bland- Altman method. 

Unlike the first two approaches which are not affected by differences in 
units or the nature of results, the Bland-Altman method analyzes the 
agreement between two estimates of the same variable [S]. Thus, we 
used SVscSA instead of INVSi, and moreover, SVscsA variations to 
compare them to SVpicco variations (A := SVVpicco — SWscSA 
with SWpicco = SVpiccoip) - SVpicco{n - 1) and SVVscsA = 
SVscsAin) — SVscSAin — 1))- The mean difference between SVVpicco 
and SVVscsA (mean A) is plotted against the average of the two volume 
variations. Mean A which represents the bias between the two methods, 
and the 95% confidence interval ([C/i„/A C/s„pA]) gives the variation 
of the values of one method compared to the other. 

3 Results 
3.1 Patients 

The 21 patients recordings were analyzed over 900 seconds, representing about 
30 averaged values, except one recording, analyzable only for the first 16 aver- 
aged values. In order to illustrate the main results, we choose the first three 
subjects in Table [TJ representing various conditions: subject one was submit- 
ted to PEEP changes (figlS]), subject two was submitted to noradrenaline dose 
changes (figlHl), subject three had no change in ventilatory parameters nor in 
drugs (figCni). 
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Table 1: Cross-correlation and linear regression coefficients between SVpicco 
and INVSi 



Subject 


Number of 
values 


Cross 
correlation 


Linear 
regression 


1 


30 


0.99 


0.98 


2 


30 


0.97 


0.97 


3 


30 


0.97 


0.97 


4 


30 


0.96 


0.96 


5 


30 


0.96 


0.96 


6 


30 


0.95 


0.95 


7 


30 


0.95 


0.95 


8 


30 


0.93 


0.94 


9 


30 


0.91 


0.91 


10 


30 


0.91 


0.89 


11 


30 


0.90 


0.89 


12 


30 


0.90 


0.90 


13 


30 


0.86 


0.85 


14 


30 


0.85 


0.77 


15 


30 


0.85 


0.86 


16 


30 


0.85 


0.85 


17 


21 


0.84 


0.87 


18 


30 


0.82 


0.83 


19 


16 


0.82 


0.82 


20 


30 


0.77 


0.77 



Tlie subjects (col. 1) are listed in the decreasing order of cross-correlation coefficients. 
Col. 2 indicates the number of measurements delivered every 30-second per subject, 
according to the PiCCO protocol and representing a 900 seconds analysis. Col. 3 
and 4 stand for coefficients of cross-correlation (mean: 0.90 ± 0.01 at lag 0) and linear 
regression (mean: 0.89 ± 0.01). Subject 21 and the last third of Subject 17, whose 
INVSi and SVpicco were not correlated at all (coefficient < 0.1), were discarded 
from the table. The last part of Subject 19 was discarded before analysis because of 
a PiCCO dysfunction. 

3.2 Cross-correlation analysis 

Table \T\ column three, shows the coefficients of cross-correlation in decreas- 
ing order. The amount of well correlated measures represents 94% of the all 
recordings. Because of an obvious divergence between SVpicco and INVSi 
(coefficients < 0.1), Subject 21 and the third part of Subject 17 were discarded 
from the table. Thus, they represent only 6% of discrepancy among all the 
recordings. Fig|3]and fig (5] respectively represent the time series of HR, MAP, 
PVR, INVSi and SVp.cco for subjects 21 and 17. 

Fig [S] represents the cross-correlation coefficients of the 20 remaining subjects 
(dashed lines). As homogeneity was verified, averaged values were also plotted 
(— • — ). The vertical axis depicts the correlation coefficients. The horizontal 
axis depicts the lag, in number of 30-second averaged values of one time series 
on another. The correlation coefficients of the unlagged data are plotted at hori- 
zontal midpoint (lag 0) . The two symmetrical continuous lines represent critical 
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Table 2: Bland- Altman test results 



Subject 


MinA 


ClinfA 


Max A 


CIsupA 


1 


-.0056 


-.0074 


.0048 


.0077 


2 


-.0026 


-.0037 


.0034 


.0039 


3 


-.0030 


-.0036 


.0022 


.0037 


4 


-.0043 


-.0057 


.0043 


.0056 


5 


-.0034 


-.0039 


.0044 


.0040 


6 


-.0048 


-.0052 


.0045 


.0053 


7 


-.0059 


-.0072 


.0043 


.0072 


8 


.0005 


-.0009 


.0009 


.0009 


9 


-.0037 


-.0051 


.0046 


.0050 


10 


-.0047 


-.0055 


.0042 
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A = SVVpicco - SVVscsA with SVVpicco = SVpiccoin) - SVpicco{n - 1) and 
SVVscsA = SVscsA{n) — SVscsAin—l). Min{A) and max{A) stand for the minimal 
and maximal values of A respectively. CI is 95% confidence interval of the agreement 
limits of A; [ClinfA CIsupA]. ★ means that all the values are inside [ClinfS ClaupA] 
except 1 value greater than CIsupA for two subjects. 



r values for a level 5 of probability (Bravais-Pearson table) . The greatest cor- 
relation stands at lag for all the remaining subjects (mean correlation= 0.90; 
sem — 0.01; p — 0.00001). This result shows an excellent temporal similarity 
between the successive measures, indicating that they change in the same way 
over time. 



3.3 Linear regression 

Table [U column four, shows the R coefficients of linear regression. Cross- 
correlation and R coefficients are strongly correlated (0.95 at Spearman rank 
order correlation test). The mean coefficient, equal to 0.89 shows a great degree 
of linearity between the two methods. FigIT] shows the plots of linear regression 
for the first three subjects. 



3.4 Bland- Altman method 

FiglT] shows the Bland- Altman plots for the first three subjects. Differences in 
the two mean SVV {SVVpicco — SVVscsa) are plotted against the mean of 
the two SVV ( Pi-CCO + SCSA_y Continuous lines respectively stand 
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Figure 4: Cardiovascular time series of subject 21, submitted to ventilatory 
and pharmacological changes. The lack of correlation between INVSi and 
SVpicco (cross-correlation coefHcient=0.039) is in favor of INVSi. Before 
increasing PEEP, nothing happens, HR, MAP, PVR are stable, thus no SV 
change can be expected. Nevertheless, SVpicco decreases then increases while 
INVSi remains quite stable. During increasing adrenaline, increasing PVR is 
accompanied, as expected, by increasing MAP and decreasing HR. An increase 
in SV is also expected, which is done by INVSi while SVpicco remains quite 
stable. 



for the mean of differences A between the two results, and the 95% confidence 
interval. Each dot stands for the difference between SVV measured by the two 
methods. In the three cases, mean A is equal to and all values are within 
the two 95% confidence intervals. Table [2] gives results for all the subjects. 
When comparing TOm(A) to C/i„/A (columns 1 and 2) and max{/S) to Clgup 
(columns 3 and 4), one can see that all values are inside [Clmf CI sup] except 
one value greater than CIsup^ in two subjects (*). These global results allow 
us to conclude that a good proximity exists between the two methods with the 
same order of dispersion and that they are interchangeable. 
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Figure 5: Cardiovascular time series for subject 17, submitted to ventilatory 
and pharmacological changes. INVSi and SVpicco are strongly correlated 
(cross-correlation coefficient = 0.84) during period A, and they arc divergent 
during period B. Decreasing noradrenaline is naturally accompanied by decreas- 
ing PVR and MAP and increasing HR. A decrease in SV is also expected, 
which is done by INVSi while SVpicco strongly increases. 



4 Discussion 

A new method for a simple and minimally invasive SVV estimation from ABP 
measurements has been validated in this study. The ABP signal is reconstructed 
with a semi-classical signal analysis method SCSA which enables the decomposi- 
tion of the signal into its systolic and diastolic parts. Some spectral parameters, 
that give relevant physiological information, are then computed and especially 
the first systolic invariant INVSi, given by the area under the estimated sys- 
tolic pressure curve. Thus, we have shown that INVSi yields reliable SVV 
assessment. 

So, in order to validate this approach, we compared INVSi estimated from 

ABP measurements with SV measured with a reference method: the PiCCO 
technique. Three statistical methods were applied for this validation: cross- 
correlation analysis, linear regression and the Bland-Altman test. Among the 
315 minutes duration of all the 21 recordings, 94 % presented a very high cor- 
relation between INVSi and SVpicco- The mean coefficient was equal to 0.9 
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Figure 6: Cross-correlation between the 30-second averaged SVpicco and 
INVSi values in 20 subjects Each line (dashed) stands for a subject; the strong 
line (— • — ) represents the average of the 20 subjects; continuous lines represent 
critical r values for a 5% level of probability (Bravais- Pearson table). The verti- 
cal axis depicts the correlation coefficients. The horizontal axis depicts the lag, 
in number of 30-second measures, of one time series on another. The correlation 
coefficients of the unlagged data are shown at horizontal midpoint (lag 0). All 
time series are exactly synchronized, with a mean correlation coefficient equal 
to 0.9 at lag (p=0.00001). 



for cross-correlation, and equal to 0.89 for linear regression. The remaining 6% 
without correlation, concerned two subjects: all of subject 21 and the last third 
of subject 17. This discrepancy, interpreted with the help of the synchronized 
HR, MAP, PVR time series, can be explained by the following remarks: 

• For subject 21 (figH]), before increasing PEEP, nothing happens, HR, 
MAP, PVR are stable, thus no SV change can be expected. Never- 
theless, SVpicco decreases then increases while INVSi remains quite 
stable. During increasing adrenaline, increasing PVR is accompanied, as 
expected, by increasing MAP and decreasing HR. An increase in SV 
is also expected, which is done by INVSi while SVpicco remains quite 
stable. 

• For the last third of subject 17 (figlH B), decreasing noradrenaline is nat- 
urally accompanied by decreasing PVR and MAP and increasing HR. A 
decrease in SV is also expected, which is done by INVSi while SVpicco 
strongly increases. 

The divergence between INVSi and SVpicco is in favor of INV Si for these 
two subjects. 

On the 94% recordings with well correlated INVSi and SVpicco, the 
Bland- Altman test showed a very good agreement between the two approaches 
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Figure 7: Linear regression plots (at the top) and Bland- Altman plots (at the 
bottom) for the three first subjects. The coefficient of correlation is greater 
than .95 for each of them, meaning a strong linear relation. All the differences 
in SVVpicco and SVVscSA are included in the 95 % confidence interval for 
each of them, meaning a good agreement between the two methods. 
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Figure 8: Cardiovascular time series of subject 1, submitted to Positive End 
Expiratory Pressure changes. A: decreasing PEEP; B: increasing PEEP; C: 
decreasing PEEP; D: increasing PEEP. INVSi and SVpicco are strongly cor- 
related (cross-correlation coefficient = 0.99). 



and demonstrated their interchangeability. It is worth noticing that this agree- 
ment is obtained in unstable hemodynamic and/or noisy conditions which prove 
the robustness of the SCSA method. Several great ventilatory or pharmacolog- 
ical changes are illustrated in figEl fig|H]and figlTU]but also by the non averaged 
time series in fig llll A noisy condition is illustrated by subject 3 on the right. 
Despite a raw PI (top right) disturbed by extra-systoles and artefacts, INVSi 
(bottom right) is well estimated. 

Therefore, this study shows that SCSA is a reliable method for SVV assess- 
ment and more suitable in the two cases of divergence. The good agreement 
between the two approaches could be explained by the fact that the main idea 
in the SCSA technique is quite similar to the PiCCO and consists in using the 
area under the systolic part of the pressure curve. However, the detection of the 
end systole with the SCSA is different from the pulse contour approach. Indeed, 
while the pulse contour approach uses an algorithm to detect the dicrotic notch, 
the SCSA uses an ABP model based on solitons that takes into account nonlin- 
ear phenomena, as described in section Il.C.l.b. This difference could explain 
the greater reliability of SCSA when discrepancies between the two approaches 
appear. Moreover, this explanation agrees with the observation of the raw ABP 
signal for subject 17: its shape is very different between the first and last part 
of the recording. 

Finally, unlike the PiCCO technique which needs periodic calibration by 
a thermodilution technique, SCSA is easier to use, requiring less equipments, 
only for ABP measurements. It is much less invasive and could be totally 
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Figure 9: Cardiovascular time series of subject 2, submitted to a vasoactive 
drug. A: increasing noradrenaline; B: decreasing noradrenaline. INVSi and 
SVpicco are strongly correlated (cross-correlation coefficient — 0.97). 




Figure 10: Cardiovascular time series of subject 3, without any change in ven- 
tilatory or pharmacological condition. INVSi and SVpicco are strongly cor- 
related (cross-correlation coefficient = 0.97). 
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Figure 11: Pulse interval (PI), Systolic Blood Pressure (SEP) and INVSi time 
series for the three first subjects with a strong correlation between SVpicco 
and INVSi . The INVSi is precisely estimated despite unstable hemodynamic 
conditions (subjects 2 and 3) and noisy conditions such as extra-systoles and 
artefacts (subject 3). 



noninvasive if we used a FINOMETER device, for instance. This latter point, 
already experimented in our previous studies [21 , should be a new perspective 
for a simple non-invasive SVV assessment. 
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